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ABSTRACT 

A solution aathod has boon davalopad for calculating coapraanlblo Invlacld 
flow through a llnaar cascade of arbitrary blada shapoa. Tha aathod uaas 
advanced surface alngularlty fornulatlons which wars adapted froa those found 
In currant external flow analyses. Tha resulting solution technique provides 
* last flexible calculation for flows through turboaachlnary blade rows. Tha 
solution aathod and aoae axagqilea of tha aathod 'a capabllltlas are prasantad. 


NONENCUrdRE 


6 

cascade stagger or blade setting angla 

As,Ba,Av,Bv, 


A 

element or panel chord length 

Ad.Bd 

source, vortax, and vortex junp Influ- 




ence coefficients 

(>n 

coordinate system associated with caa- 




case elements 

C 

blade chord 

o.y.4 

source, vortex, and vortex Jump alngu- 

c 

curvature paraaeter 


larlty strengths 

T.j 

unit vectors 

* 

total velocity potential 

t 

a control point on the cascade surface 

d 

component of total velocity potential 

n.I 

surface unit nornal and tangential 

e 

orientation angle of panel or velocity 


vectors 


point coordinate aystam with respect 
to reference coordinate system 

P 

arbitrary point In apace 

p 

fluid density 

S 

cascade pitch or gap . 

V 

gradient vector operator 

8 

surface distance from center of panel 

Subscripts: 


V 

Velocity relative to the blade 

l.c 

Incompressible and compressible 

V 

disturbance velocity 

E.I 

exit and Inlet conditions 

X.Y 

r(x - O/S, »(y - n)/S 

l.j.k 

panel Indices for source, vortex, and 

x.y 

coordinate system associated with the 


vortex Jump alngularlty strengths 


point trhere the velocity Is to be 
calculated 

N,T 

normal and tangantlal component on 




blade surfaces 

a 

flow angle relative to axial direction 

TE 

trailing edge quantity 



t 

stagnation or total condition 



U,L 

upper and lower surfaces 



onset 

mean undlstrubsd flow 


Ovarbara: 


( ) 

vector quantity 

mi 

( ) 

averaged quantity 


P ♦ f PV* 


conatant ■ P^ 
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IMTROmjCTION 

Currant daalgn of turboaachlnary blada rowa ra- 
lia? n tha uaa of coaputar codaa to nodal tha flow 
or H. da-to-blada aurfacaa. Moat of tha codas cur- 
rv -*'/ uaad for blada daalgna nodal tha flow aa 
Invlacldi Irrotatlonal and conprcaalbla. Tha flow 
flald la solved using flntta dlffaranca or flnlta 
alanant nuasrical tachnlquaa. Tha aoluttona glvan 
by thasa codas can bo quite accurate. However, It 
has been the experlenca of the author that use of tha 
field techniques can require an experienced user to 
SMnlpulate Input data and control paraaetara. and 
that these codes nay also Halt a designer In the 
types of blade geonetrlea, cascade configurations, 
and flow condltlona that can be considered. The 
designer can therefore spend a considerable portion 
of hla effort In aaklng a prellnlnary design confora 
to the raqulraaents and eccentricities of a particu- 
lar flow analysis code. In order to accelerate the 
blade design process and give the designer aore free- 
doa In developing blade shapes, s slnple blade-to- 
blade flow code has been developed. The objective of 
the code Is to provide a versatile, stable, and effi- 
cient calculation which has sufficient accuracy to 
detaralne the basic aerlts of a particular blade row 
design. The slnple coda provides a means of rapidly 
screening blade designs for analysis by nore cooiplex 
codes which have better accuracy. 

The panel or Integral equation solution tech- 
nique was selected as a baals for the slnple design 
code. This technique neets the requirements of being 
a versatile, stable, and efficient calculation 
acheaw. It is however an Incompressible solution and 
requires the use of a compressibility correction to 
account for compressible flow effects. 

Panel solutions have been used for several years 
by external aerodynamlclsts who have developed and 
refined them Into one of the primary design tools of 
the aircraft industry. In the past. Internal flow 
aerodynamlclsts such as Hartensen [l]. Gelsing [ 2 ], 
and Hess and Stockman [l] have made use of Integral 
solutions, but the basic technique has not been pur- 
sued extensively for Internal flow problems. 

In this paper a panel solution Is described for 
Inviscld cooqiresslble flow through a plane cascade. 

An Incompressible solution is developed using 
advanced Integral equation formulations. Compress- 
ibility effects and stream channel varlatons are 
Included by adapting the Llebleln-Stockman correc- 
tion technique ^rom engine Inlet calculations 
to the cascade problem. Several example aolutlons 
are given to demonstrate the capabilities of the 
method. 


ANALYSIS 

The paneral cascade problem Is shown In f Ig- 
ure 1. The governing equations and boundary condi- 
tions for an Invlacid Incompreaslble caacade flow are 
the following: 

V • V - 0 (1) 


*^body 

eurface 

autfece 

(3) 


V Vj aa 

X 

(4) 


V Vg aa 

X 

(5) 


The solution Is developed using a velocity potential 
which Is the sum of a conatant velocity potential plua 
a disturbance. The conatant velocity potential la 
determined from the swan of the upatreom and downstream 
velocities. The disturbance quantities ace unknown. 
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The constant or onset velocity and potential satisfy 
continuity (eq. (1)) exactly. The disturbance poten- 
tial and velocity are to be determined so that contin- 
uity Is satisfied. 


V • v " V ^4 - 0 

(7) 

^onset ■ 

(8) 


The flow field la determined by solving equations (6), 
(7). and (8) subject to boundary conditions (3), (4), 
and (5). The laonantum equation (2) does not enter 
Into the calculation except to relate the pressure to 
the velocity field. 

The equation for the disturbance potential (eq. 
(7)) Is Laplace's equation. Since It is a linear 
equation, superposition of known simple solutions of 
Laplace's equation may be used to develop nore complex 
solutions. A general solution to flow over a body or 
cascade of bodies may be developed by using basic in- 
compressible potential flow solutions for source and 
vortex flows distributed along the body surfaces, and 
varying the strength of the source and vortex singular- 
ities so that the problem's boundary conditions are 
satis! led. This Is the basis of the panel method as 
described by Hess and Smith [sj. 

In practice the surface of the body Is represented 
by inscribing a polygon as shown In figure 2. The 
simple source and vortex singularity solutions are de- 
scribed over each element, and a control point Is 
selected on each element where .he normal velocity 
boundary condition Is to be applied. There will be n 
element end points and n - 1 control points. The 
first and nl^ element end points are coincident. 

The simplest form of the solution would be to use 
flat elementa with constant source and vortex slnular- 
ity strengths over each element and the control points 
located at the center of the flat element. However, 
Hess [b] demonstrated a means of developing higher 
order accurate panel solutions by using curved elements 
with varying singularity strengths and locating the 
control points closer to the true body surface. Build- 
ing on Hess's work, a higher order accurate panel solu- 
tion Is developed for linear cascades. 
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C «>c«d> p»n«l forauUtlon . A alngl* caacad* 
panal la ahown in flgura 3. Tha panal haa aaaoci- 
atad with_U a cartealan coordtnata ayataa with unit 
vactora 1,J. Tha eoordlnata ayataa la cantarad on 
tha panal. Ita origin la loeatad at tha panal con- 
trol point and tha x or ( axla la aligned with 
tha panal chord. Tha velocity at any point P due 
to tha alngularlty diatrlhutlon along each caacada 
panal la to ba dataralnad. Tha coordinates for 
point P are denoted aa x.y and thoae for polnta 
on tha panal aa C,n. 

Heaa ■ bowed that the velocity at any point P 
due to a alngle panel with a aource atrength diatrl- 
butlon, o(a)i la given by 

i/2 

V - 2 J I + J^o(a) || dC (9) 

-A/2 ^ 


whera 


r • y<x - * (y - n)^ 

For a caacade problco a aeries of these integrals is 
needed. 
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where 

rf • ^(x - O* + y* 


Tha factor c appearing in tha aeriaa expanalon ia 
related to tha panal curvature and ia deteminad by a 
throe point curve fit uaing the averaged alanant alopaa 
at the end pointa of the panel. The expanaion of 1/r’ 
ia needed for nuanrical atability. Aa point P ap- 
proaches tha path of integration tha integrand baeonaa 
very large and the reaulting axpreaaion for velocity 
bacoaias divergent. The expanalon in tama of rf 
alleviataa thla problen. 

In the praaent atudy equation (11) waa integrated 
nunerlcally. The aarlea axpanslona for q«aiO in 
terns of i were substituted Into the integrand. Tha 
resulting expression was evaluated using a three point 
(C • - A/2.0.A/2) Slnpson's rule. Again, nunarical 
difficulties %nre encountered In calculating tha inte- 
grands as the velocity point approached the path of 
Integration. However, fron aquation (10^ it can ba 
aeen that only ths first tern In tha serif's Is causing 
the difficulty and that this tern la Identical to the 
velocity equation for an Isolated body (eq. (9)). 
Therefore, the nunerical problen nay be avoided by sub- 
tracting the first tern of the series fron tha inte- 
grand, Integrating, and addlug Hass's analytical 
expression to the Integral. Equation (11) becomes 
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slnh(X)cosh(X)l * sln(Y)cos(Y)1 


slnh^(X) + sln^(V) 
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(X - Qii-(y_- Cn - hs))I 
(X - 0^+(y - (n - ks))2 



(oo + oiCHDdc + 


^Isolated 

body 


(13a) 




(X - c)i » (y - ( n kS ))j 
(X - £)2 + (y - (n + kS))2 
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where k Is the sunuation iedex for cascade bodies. 
This aeries Is the caacade potential velocity which 
was shown in reference 2 to be 
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v-li 
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where 


It 

-A/2 


slnh(X)cosh(X)t ■»eln(Y)cos(Y).1 


slnh^(X) + sln^(Y) 


^o(s)||de 


( 11 ) 


X - I (X - t); Y - I (y - n) 


Hess ^6^ developed the higher order panel for- 
mulas by Integrating equation (9) using series expan- 
sions for the different terms in the Integrand. 

These expansions for a second order accurate panel 
solution are 

Ti - c£^ + oah 

O - Oq + ois ♦ O(s^) 

a - £ + I c^(2 + o(£*) 

-L . -I 


ff 
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Equation (13a) has been developed for a source 
singularity distribution. It can also be used for the 
vortex singularity distribution since from incompress- 
ible potential flow theory the source and vortex sing- 
ularities are related. The vortex equation would be 
similar to equation (13a) except o w>uld be replaced 
by Y> the velocity components change. The resulting 
equation Is 


uY)cos(Y)l - slnh(X)cosh(X)l 


slnh^(X) + sln^(Y) 


■ {W)]' 


(Vo + YiUa)d£ 


(13b) 


Integral equation fomulatlon . The Integral equa- 
tion formulation of Bristow [7J was used In the present 
cascade panel solution. The primary feature of this 
formulation Is the use of Green's third Identity [sj. 
Bristow denonstrsted that by using this identity for 
two dimensional flow a direct relationship can be de- 
veloped between the onset flow, and the source and 
vortex singularity distributions. Along the body sur- 
face the tangential and normal velocity will be given 
by 
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Tb«a« ralatloaahlpa Man Chat tha ala|ularlcy 
atraagtha will raaatn boundad ragardlnaa of tha flow 
boundary condttlona or gooMtry. Brlatow [71 ahowad 
Chat tha nlld atngularlty atrangth dtatrlbutlon gan- 
aratad by ualng Graan'a third Idantlty will raault 
In panal aolutlona which ara laaa aanaltlva to coor- 
dlnato awoothnaoa and thin airfoil gaoMtrlaa. 

Baforo awklng uaa of Graan'a third ldantlcy> an 
approalaaclon of aquation (6) la naadad. Thla can 
bo wrltcan In acalar conponanta aa 
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The first tarns In the equations (IS) are the 
contribution of the constant or onset velocity to the 
total velocity vector in equation (6). The angles 6 
and a are the orientations of the coordinate system 
associated with point P and the onset velocity vec- 
tor with respect to a reference direction. The 
suanatlon terns represent the disturbance velocity 
V In equation (6). Each tens in the aunnatlon is 
an Integration of equations (13a) or (13b) over one 
of the elenents which sake up the cascade body. The 
A's and B's in the sunnatlons are coefficients of 
the singularity strength values fron the elenent 
Integrations. These Influence coefficients depend on 
the geometry of the problem and must be recalculated 
for each point where a velocity calculation is to he 
made. The source strengths a which are assumed to 
vary linearly over the element surface are defined 
In terms of the element control points using a three 
point curve fit. The vortex atrengtha y also are 
assumed to vary linearly, but they are defined In 
terms of the eleewnt end points using four point 
divided differences. The additional singularity 
strengths S are part of the vortex singularity 
distribution. They represent a Jump in the vortex 
strength which Green's identity, equation (14a), 
Indicates will occur st discontinuities in surface 
slopes. Such discontinuities are encountered for 
example at sharp leading and trailing edges of 
blades. When a surface slope discontinuity Is en- 
countered In the calculation, ena sided curve fits 
and differences are used for the source, vortex, and 
geometry terms on each side of the discontinuity. 

Normalising equations (14) and (IS) by the onset 
velocity msgnltude, applying equation (13) to a con- 
trol point 1 on the body surface, and resolving tha 
calculated velocity Into tangential and normal veloc- 
ity components, expressions for tha surface velocity 
nay be written which relata Green's third Identity 
to the Integral equation approxlnatlon, 


■ - ala(e| - n) ♦ 


(16b) 


B ft-l ■ 

V- - com(B^ - a) + V AsgiOi ♦ T ♦ V 

1-1 J-1 k-1 (16c) 

n-1 

Vw 

» - Pi " ■ ”l " Pi 

(16d) 


n n -1 B 

Vq- - •*n(6i - a) ♦ 5) Bagioi + Bvtjtj ♦ ^ Bdik<t 
* l5i J-1 k-1 


where y with an * denotes that tha vortex strength 
la evaluated at a penal control point while y with- 
out an * la evaluated at the penal and point, 

A ayatam of llnaar algebraic equatlona to calcu- 
late the singularity atrengtha may now be constructed. 
This system of equatlona la sinllar to those developed 
by Bristow [?3 for Isolated bodies. These aquations 
are generated by using tha noimal velocity surface 
boundary condition, farstream boundary conditions, and 
a tangential velocity error mlnlmlaatlon aquation. 

The normal velocity boundary condition la first used 
with equation (16b) to solve for the source strengths, 
and then with equatloo U6d) to give n-1 equa- 
tions. Tha farstream boundary conditions, tha given 
inlet and exit flow angles, are used to calculate the 
circulation of the cascade. This la used with equa- 
tion (16c) to provide one equation. The tangential 
velocity error minimisation equations are needed if 
there are surface slope discontinuities, and vortex 
Jump strengths are to be calculated. One minimisa- 
tion equation will be added to the system for each 
surface slope discontinuity. The tangential velocity 
error function is the sum of the squares of the dif- 
ference between equations (16a) and (16c) multiplied 
by the local panel chord length. 

This system of aquations satisfies all the prob- 
lem boundary conditions and Includes all the unknown 
singularity strengths. However, there Is one more 
equation than unknowns. This occurs becaus'. the 
first and last element end points are coincident for 
a closed body and therefore yi and Yq are equal. 
The additional unknown used to complete the system 
is an assumed uniform error In the normal velocity 
boundary conditions. The uniform error verlable is 
denoted as c and is Included in the normal velocity 
equations. 

In final form the equations to be solved for the 
singularity strengths and uniform error are: 

1. Normal velocity surface boundary condition 

Oj - sln(6j^ - o) + Vjj^ (17a) 

n m n-1 

X + ^ Wtk«k“V„ ♦•*“<6t-o) ®*tl®l 

J-1 k-1 1-1 


2, Circulation 

r - S cos(a) [tan(oj) - tan(og)^ (17b) 


+ 0 V-ds 
•'body 


(17c) 


4 



This lystM of •quatlona la aolvad ualag a dlraet 
Mthod. Tha tachniqua uaad ia known aa Cholaaky'a 
■athod. It raaolvaa tha coafflctant natrix Into uppar 
and lowar triangular natrlcaa ao that additional aolu- 
tlona for varlatlona la aurfaco nortMl valoeltlaa and 
faratraan boundary condltlona raqulra only a back aub- 
fl7d) atltutlon. Thla proparty la uaad la tha Itaratlon 

' procadura for tha Rutta condition and In tha calcula- 

tion of oHiltlpla Inlat and axlt flow angla caaaa. 


3. Tangantlal error nlnlnliatlon 



(17e) 


If only the upstrean flow conditlont are given, 
the circulation cannot be calculated directly and 
equation (17d) la replaced by a Kutta condition. An 
iteration procedure ia then used to find the proper 
onset flow direction and circulation. The Kutta con- 
dition for sharp trailing edgea is the sane as 
Bristow Qj used, that is, flow normal to the trail- 
ing edge bisector is zero. For the round trailing 
edge found in turbomachinery, this condition is mod- 
ified. Two-line segments are constructed parallel to 
the trailing edge bisector from the trailing edge 
circle tangency points of the upper and lower blade 
surfaces. At a point on each line Just downstream of 
the trailing edge, the normal flows to these lines 
are required to be equal and opposite. The two types 
of Kutta conditions are illustrated in figure 4. The 
Kutta condition equation Is written as 

n 

^ (®''TE,U ®''TE,l)j7j + 2c " I sln(e^£ - a) 

J-1 


n-1 

- ^ <®®TE.U ®®TE,L^j^''i 

1-1 

In the iteration procedure equations (17a), (17b), 
(18), and (17e) are solved using a guessed value for 
o. The resulting singularity strengths are used in 
equation (17d) to determine the circulation. The 
calculated circulation and onset flow angle a are 
then used in the following relation to determine the 
upstream flow condition 

. tan-1 ninioL;^^ 

I L cos (a) J 

If the calculated value of Oj differs from the 
specified inlet value, the onset flow angle a is 
updated and the procedure is repeated. 


Cowpresalblllty corrmctton. Tha Intarnal flow 
coapresslblllty correction of Llabloln and Stockaan 
la used in this study. This correction has bean used 
succaaafully for a nunber of years to calculate angina 
Inlet flows. The correction la axprassad In tha fol- 
lowing form 

Vi/i?i 

• Vj(pj/p^) (20) 

where la the corrected compraaslbla valoclty, Vi 
is the local Inconpresslbla velocity froa the panel 
solution, pt la the Incoaprasslble density taken to 
be the stagnation densltjr Pt.'^c the average coa- 
presslble density, and Is the average Incoapreaa- 
Ible velocity across the flow paaaage. Tha average 
compressible density is found by solving the following 
expression derived from the one-dlaenalonal continuity 
equation and the laentroplc relations. 



where y 1* the ratio of specific heats and la 

the critical velocity. The average incompressible 
velocity is calculated on planes which are nonssl to 
the mean flow and which contain the local velocity to 
be corrected. 

For inlets with little streamwlse turning of the 
flow, the application of the compreasibillty correction 
is straightforward, the mean flow direction can be 
assumed to be approximately normal to axial statlona 
and the average incompressible velocities are calcu- 
lated using the axial station flow areas. However, for 
cascade flows with significant asKnints of flow turning, 
the mean flow direction cannot be assumed and must be 
calculated. In this study, the components of the mean 
velocities are fiist calculated at axial atatlons by 
using continuity and circulation equatlona. Next the 
mean flow direction through the blade paasage is de- 
termined from the axial station velocity components. 
Planes normal to the mean flow are then constructed 
for each surface velocity to be corrected. Finally, 
the average Incompressible velocity on the normal 
planes is found by interpolation of the axial station 
values. 

Stream channel variation corrections . A simple 
correction for given stream channel variations la made 
during the compressibility correction calculation. The 
incompressible surface velocities and isean flow are 
calculated assuming a uniform stream channel height 
using the inlet flow conditions as a reference. These 
velocities are then scaled by the ratio of the inlet 
to local specified stream channel height. The reault- 
ing scaled velocities are used in the compressible 
corrections of equation (20). 
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SAMPLE CALCUtATIOHS AND DISCUSSION 

Four MBplo calculatloM ara praaantad to 
daBoaatrata tha aolutton aethod. Tabla I ahowa tha 
ganaral flow and aolutlon conditlona for aach calcu- 
lation. All anaaplaa wara calculatad ualng tha 
IBM 3033 coaputar at Lawla Raaaarch Cantar. A 
FORTRAN coaputar coda of tha aolutlon aathod waa 
davalopad for tha calculatlona. Moat blada ahapaa 
«iara approxlaatad by ualng AO to 60 alaaanta and tha 
aolutlon raqulrad laaa than A aaconda of coaputar 
tlaa. Nhlla thla coaputatlon tlaa la vary rapid. It 
ahould ba notad that coapraaalbla flald aathoda hava 
baan rwrtad which hava run tlnaa cn tha aaaa 
ordar [sl. In tha Goatalow and Sana caaaa, tha caa- 
cada airfoil ahapaa wara daacrlbad dlractly ualng 
dlacrata polnta. Howaver, In t^a Kataanla and Hobaon 
caaaa, tha caacade boJy ahapaa wara ganaratad by 
ualng apllne flta for tha uppar and lowar av.'facaa, 
and clrclas or alllpaea for the rounded laadtAg and 
trailing adgas. 

Incoaipreaslble caacade flow. The exact incon- 
preaalble solution of Goatalow Flflf] la conparad trlth 
tha present solution In figure S'. The large velocity 
gradients In the leading edge region and tha flow 
past the thin trailing edge are well Modeled by tha 
panel solution. Tha overall agreement between the 
two solutions Is very good. 

Compressible cascade flow . In figure 6 a com- 
parison is shown between the panel solution and the 
finite difference solution of Katsanls [ll^ for a 
turbine stator. Blade surface velocities are In the 
high subsonic range. The overall solution agreement 
is good, but differences do occur along the pressure 
surface and near the suction surface trailing edge. 
These differences are attrUutsble to the compraa- 
slblllty correction, since the two solutions were 
found to be Identical for Incompressible flow through 
the same cascade. 

The general accuracy of the panel method Is also 
ohown in figure 6. Two panel solutions are shown 
with a different number of panels used in each case. 
The less accurate solution with fe«fer panels differs 
slightly from the more accurate solution. Most of 
this difference occurs In the leading edge region 
where large variations in surface curvature occur. 

The coarser paneling Is unable to resolve the rapid 
change In curvature and as a result the curvature 
effect on the flow Is smeared over a larger surface 
area. By considering the terms used in developing 
the cascade panel solution, the accuracy of Che solu- 
tion should be on the order of the panel size 
squared; although this accuracy has not been 
verified. 

Transonic cascade flo w. Comparisons with tran- 
sonic hodograph solutions are shown In figures 7 
and 8. In figure 7 the comparison Is with Hobson 
solution C^2j, and In figure 8 the comparison Is with 
a Sam cascade blade design using the hodograph solu- 
tion technique of Bauer, Garabedlan, and Kom L^3j. 
The panel solution compares well with the hodograph 
solutions for the sophistication of the computation 
Involved . 

Deacrspancles In the solutions are caused pri- 
marily by the compressibility correction which cannot 
be expected to exactly pr-dlct the flow at such a 
high level rf compressibility. The large velocity 
spike at the leading edge of the compressor stator is 


due CO tha lack of blada gaomatry Information given In 
the region, and tha axtranaoua v^oeltlas downatream 
of tha blade trailing edge ara caused by an extension 
the blada downatream to account for the blade waka. 


KDTTA GONDiTItW AND STREAM CHANNEL VARIATION 

In the aanpla calculations Just presented, tha 
downstream flow angle was given. These calculatlona 
have also been made ualng the Inlet condition and a 
Kutta condition to datamlna the domatraam flow con- 
dition. The calculated downstream flow angles ware 
31.190 (Goatalow), -67.88° (Katsanls, 32). -67.23° 
(Katsanls, 98), -AA.16° (Robson), and -31.0° (Sans). 

In tha first two casas, Goatalow and Katsanls, the 
flow leaves tha blade surface smoothly with little 
deviation from the trailing edge blade angle and tha 
Kutta condition works well* However, In the Hobson 
and Sans cases there Is significant flow deviation at 
the blade trailing edge and the Kutta condition re- 
sults are poor. Soas laprovsawnte In these results 
could be niule by Including s deviation model In the 
Kutta condition formulation. 

The effect of etroam channel variation on the 
flow has not been studied extensively. Some compari- 
sons have been made trlth the finite difference code of 
Katsanls Ql3 for subsonic flows and the agreement was 
good. 


OONaUSION 

Current external aerodynamic Integral equation 
techniques uy be adopted for use In Internal flow 
calculations. The Inherent computational speed and 
flexibility of Integral equation solution can make 
them very useful for design calculations. The present 
cascade aethod Is a versatile design tool that will 
allow a designer to explore many preliminary blade 
designs In a short period of time. Although the 
BMthod does not give exact solutions for compressible 
flows, example calculations do demonstrate that it Is 
sufficiently accurate to provide a means of selecting 
blade designs for further analysis. 
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TABLE I. - SAMPLE CALCULATIONS 


Case 

Inlet 

Exit 

Cap 

Chord 

Nunber 

of 

panels 

CPU sec 
IBM TSS/370 


V/Vcr 

“I 

v/v 

“E 

Gostelow 

— 

53.5 

— 

30.025 

1.239 

49 

2.42 

Katsanls 

0.231 

0 

0.727 

-67.00 

0.747 

32/98 

1.54/7.30 

Hobson 

0.610 

46.123 

0.610 

-46.123 

0.526 

61 

3.43 

Sanz 

0.826 

44.61 

0.529 

3.91 

-1 1 

0.863 

51 

1 

2.61 
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Figure (k - Compressible cascade flow. 
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Figure a. - SugercrttkaJ compressor strtor. 



